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Energetic particle production in supernova remnants (SNRs) through diffusive shock 
acceleration has long been suggested as a mechanism by which the bulk of cosmic rays 
with energies less than 10^^ eV are formed. If highly efficient, this process can 
have a significant effect on the X-ray emission from SNRs as well as their dynamical 
evolution. Here we investigate the expected modification to the thermal X-ray emission 
from the forward shock region of young SNRs in which efficient particle acceleration 
is occurring. Using hydrodynamical simulations for a range of ambient density and 
magnetic field values, we produce spectra for both the thermal and nonthermal emission 
components of the postshock gas. For a given ambient density and explosion energy, 
we find that the position of the forward shock at a given age is a strong function of the 
acceleration efficiency, providing a strong signature of cosmic ray production. Using an 
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approximate treatment for the ionization state of the plasma, appropriate for the range 
of models considered, we investigate the effects of slow vs. rapid heating of the postshock 
electrons on the ratio of thermal to nonthermal X-ray emission at the forward shock. 
We also investigate the effects of magnetic field strength on the observed spectrum 
for efficient cosmic ray acceleration. The primary effect of a large field, aside from 
an overall increased flux for higher fields, is a considerable flattening of the nonthermal 
spectrum in the soft X-ray band. Spectral index measurements from X-ray observations 
may thus be used as indicators of the postshock magnetic fleld strength. The predicted 
gamma-ray flux from inverse-Compton (IC) scattering and neutral pion (vr*^) decay is 
strongly affected by the ambient conditions. For the particular parameters used in 
our examples, the IC emission at ~ 1 TeV exceeds that from vr*^ decay, although 
at lower energies of several GeV and at higher energies >TeV, this trend is reversed 
for cases of high ambient density. More importantly, for high magnetic fields, we find 
that radiation losses, combined with evolutionary effects, produce a steepening of the 
electron spectrum over a wide energy range and results in a severe flattening of the IC 
spectral shape, making it more difficult to differentiate from that due to pion decay. 

Subject headings: Supernova Remnants, cosmic rays, shock acceleration. X-ray emission, 
MHD turbulence 



1. INTRODUCTION 

The production of relativistic electrons is beyond doubt in young supernova remnants (SNRs) 
and it is believed that remnants simultaneously produce relativistic ions, i.e., cosmic rays (CRs), 
although the evidence for this is less direct. SNRs also produce thermal X-ray line emission, and 
this contains a vast amount of information about the composition and ionization state of the gas 
that is absent in continuum emission. The premise of this paper is that shocks in SNRs accelerate 
ions efficiently via diffusive shock acceleration (DSA) and that this particle acceleration infiuences 
the evolution of the remnant, and in particular, the thermal X-ray emission, in important and 
predicable ways. In the current era of high resolution, high energy astronomy, it has become 
critical that self-consistent, broad-band models of SNR radiation under the assumption of efficient 
particle acceleration be developed. Towards this goal, we present a self-consistent model of thermal 
and non-thermal X-ray emission in SNRs. The model includes the acceleration and feedback of 
CRs on the SNR evolution and on the structure and thermal properties of the shocked gas in the 
interaction region between the forward and reverse shocks. In this paper, the first in a series, we 
investigate the dynamical and spectral evolution of the forward shock region only, where the vast 
majority of the cosmic ray production is expected to occur. In subsequent papers we will treat 
results for the ejecta component and provide a more complete treatment of the ionization within 
the SNR. 
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Radiation from shell-like SNRs consists of thermal emission from shock-heated gas and non- 
thermal emission from shock accelerated particles^ X-rays consist of a thermal component which 
contains emission from both shock-heated ejecta and circumstellar medium (CSM). Since the ejecta 
and CSM generally have very different compositions, the separation of these two components 
observationally is often possible. In this paper we concentrate on thermal emission from the shocked 
circumstellar medium (CSM) material in order to highlight the physical effects we wish to emphasize 
without the additional complications introduced by the complex composition and ionization state 
of the shocked ejecta. 

A number of SNRs also show non-thermal X-ray emission and some (e.g., RX J1713. 7-3946) 
even show non-thermal X-rays with no discernable thermal component (|Slane et al.lll999l ). Non- 
thermal X-rays are believed to be synchrotron radiation from shock accelerated TeV electrons. 
Radio emission has long been observed from SNRs in the form of synchrotron radiation from 
relativistic electrons accelerated at the outer blast wave. In addition, higher energy photons 
at GeV-TeV energies (either inverse-Compton radiation from electrons or pion-decay emission 
from ions) have been detected from some re mnants, particularly by HESS and CANGA ROO (e.g., 



Aharonian et al 



2004; 



VERITAS (ILeBohec et al.l 120061 ). 



Tanimori et al.lll998l ) and now MAGIC (e.g., I Albert et al.l 12009 ) . and soon 



If the particle acceleration process is efficient DSA, as is now becoming more generally accepted, 
the nonlinear theory makes clear predictions for the underlying electron and ion spec tra which 



produce the broad-band photon emission (e.g.. lJones &: Ellison 



1991 



Baring et al.lll999l ). While a 



number of authors have calcula ted the continuum emission from synchrotron, inverse-Compton 



2001 



Berezhko. Ksenofontov. fc Volk 



and pion-decay emission (e.g., lEllison. Slane. h Gaensler 
2002I ). there have been, to our knowledge, no self-consistent calculations of the thermal and non- 
thermal X -ray emission in SNRs. The most co mplete attempt to do this that we are aware of 
is that of iDecourchelle. Ellison, fc Balletl (j2000l ). who used a self-similar analytic model of the 
SNR hydrodynamics with input from a nonlinear calculation of DSA, which modified the thermal 
properties of the interaction region. This model did not, however, self-consistently calculate the 
S NR evolution or inc l ude th e non-th ermal synchrotron component self-consistently. Earl ier work 
by lDorfi fc Bohringen (|l993l ) (see also lChevalieiill983l : lHeavenslll984l : lBoulares fc Coxlll988l ). clearly 
showed the importance of including nonlinear particle acceleration in SNR models but did not 
attempt self-consistent calculations of thermal and non-thermal emission in the X-ray band. 

In this paper, we have incorporated a non-equilibrium ionization (NEI) model of X-ray line 
and continuum emission into a spherically symmetric hydrodynamic simulation of SNRs. The 
hydrodynamics include efficient, nonlinear DSA and the effects of this particle acceleration on 
the shocked thermal plasma, as well as on the structure and evolution of the remnant. Our 
model includes: (i) shock heating at both the forward and reverse shocks; (ii) nonlinear particle 
acceleration at the forward shock; (iii) a consistent modeling of the evolution of the material in the 



^We do not consider emission from heated dust or from a compact object in this paper. 
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interaction region between the forward and reverse shocks; and (iv) the calculation of full electron 
and proton distribution functions which allow the non-thermal synchrotron contribution to the 
X-ray emission to be determined self-consistently with the thermal emission. The thermal emission 
is calculated for an optically thin plasma in nonequilibrium ionization with negligible cooling. 

Efficient DSA predicts greater compression ratios and lower plasma temperatures than produced 
in shocks with little or no particle acceleration. These two effects will influence the X-ray emission 
from the shocked CSM where DSA is expected to be efficient. The changes in SNR evolution 
produced by efficient DSA at the forward shock (FS) modify the thermal emission from the shocked 
ejecta even if DSA is absent at the reverse shock (RS). While these changes may be small, the 
non-thermal synchrotron emission from the FS may produce large changes in the observed ejecta 
emission through line-of-sight effects even without DSA at the RS. 

In § 2 we outline our model. We use the results from time-dependent hydrodynamical models 
with efficient DSA as inputs to a NEI calculation. We discuss our approach to the problem of 
nonequilibrium ionization in § 2 as well. We present our results in § 3 and discuss the quantitative 
and qualitative effects of DSA on thermal X-ray emission in § 4. Finally, in § 4 we present and 
discuss the limitations of our approach and present possible solutions to these limitations as well 
as future directions. 



2. CR-HYDRO MODEL WITH X-RAY EMISSION 

Our model consists of three main parts: (i) the hydrodynamics of the SNR coupled to nonlinear 
DSA (our so-called CR- hydro model); (ii) a calculation of X-ray synchrotron emission from shock 
accelerated electrons; and (iii) a non-equilibrium ionization (NEI) calculation of thermal X-ray 
emission using the output of the CR-hydro model. 



2.1. CR Hydrodynamics and Electron Spectra 



We calc ulate the hydrodynamic evolution of a SNR with a radially symmetric model described 
in detail in lEllison &: Cassam-Chenaj (120051 ) and references therein. For reference, we reproduce 
the geometry of the simulation in Figure [H The current model couples effici ent DSA to the 
hydrodynamics and differs from that des cribed in lEllison &: Cassain-Chenail (120051 ) in that we have 
replaced the alge braic model of DSA of Iserezhko fc ElhsonI dlQQfll) with the more accurate semi- 
analytic model of Blasi. Gabici. &: Vannoni ( 20051 )1^1 Given an injection parameter, Xinj (this is ^ 
in equation (25) in 



Blasi et al 



200 



B), the semi-analytic model solves the nonlinear DSA problem 



The work of 



Blasi et al. I (|2005l ). which is based on the approximate solution of iBlasil (|2002l . |2004| '|. showed that 



multiple solutions could be eliminated by relating the injection momentum to the th ermal momentum of the shoc ked 
particles. Recently, an exact solution for arbitrary conditions has been presented by Amato fc Blasil ( 20051 . 2006h . 
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Fig. 1. — Spherical SNR model showing the forward shock (FS), reverse shock (RS), and contact 
discontinuity (CD). The CD separates the shocked circumstellar medium (CSM) from the shocked 
ejecta. The RS is an inward directed shock even though it moves outward during the early stages 
of the SNR evolution. As the forward and reverse shocks overtake new material, the cosmic-ray 
population is determined in concentric shells (indicated schematically) that evolve as the SNR ages. 
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Table 1: Parameters for models shown in Figures [2H5] and Figures [7Ml2[ In all models Mcj = 1.4Mq, 
*SNR = 500 yr, -Dsnr = Ikpc, and an exponential ejecta density profile is used. For all models 
except D2, E^sN = 10^^ erg. For D2, E^sN = 1-4 x 10^^ erg. 



ModeP 


nn 


Bo 


Xinj 


ercl 


<^csc 


not" 


B2^ 




/ synch \f 
\ therm ' 


pion 1 

T^llTeV 




[cm~^] 






[%] 


[%] 






[pc] 






Al 


0.1 


15 


6 








4 


60 


5.0 







A2 


0.1 


15 


3.91 


10 


1.6 


4.3 


65 


4.9 


4.68 


0.048 


A3 


0.1 


15 


3.62 


50 


16 


6.5 


98 


4.7 


13.8 


0.20 


Bl 


0.01 


15 


3.58 


36 


7.5 


5.5 


83 


6.6 


3900 


0.010 


B2g 


0.1 


15 


3.73 


36 


9.4 


5.5 


83 


4.8 


13.5 


0.077 


B3 


1 


15 


3.82 


36 


11 


5.6 


84 


3.4 


0.44 


0.65 


CI 


0.1 


3 


3.752 


36 


11 


5.6 


17 


4.9 


1.5 


0.027 


C2g 


0.1 


15 


3.73 


36 


9.4 


5.5 


83 


4.8 


13.5 


0.077 


C3 


0.1 


60 


3.2 


36 


19 


5.2 


310 


4.7 


22.7 


1.1 


Dl 


0.1 


15 


3.3 


63 


23 


7.7 


116 


4.6 






D2'' 


0.1 


15 


3.3 


64 


24 


7.8 


120 


5.0 







"Models A, B, and C are shown in Figures [S] |9l and 1101 respectively. Models Dl and D2 are shown in Figure [2] 
''This is the percent of energy flux which escapes upstream of the FS at tsNR. 
''This is the FS compression ratio at tsNR- 

''This is the magnetic field immediately behind the FS at tsNR- 
■^Radius of the FS at tsNR- 

■'^These values are summed over the 0.4 — 10 keV energy range and are for the Coulomb heating case. 
"Models B2 and C2 are identical. 

'^This model has Esn ~ 1-4 x 10^^ erg and should be compared to model Al with Esn = 1 x 10^^ erg. 
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at each time-step of the hydro simulation using the shock speed, shock radius, ambient density 
and temperature, and ambient magnetic field determined in the simulation. With the accelerated 
particle distribution, an effective ratio of specific heats is calculated and used in the hydrodyii amic 
equations, completing the coupling between the two (see lEllison. Decourchelle. &: Balletl |2004| . for 
a fuller discussion). 

At any given time, the CR-hydro calculation returns a set of concentric shells of shocked 
material, both shocked ejecta and shocked CSM, separated by a contact discontinuity (CD). The 
hydrodynamics are calculated in Lagrangian mass coordinates so the amount of material in a mass 
shell remains constant as the volume of a shell is adjusted in response to the changing pressure. 
New shells of shocked gas are added as the forward and reverse shocks overtake fresh material. 

Particle acceleration influences the SNR evolution since relativistic particles produce less 
pressure for a given energy density than do non-relativistic particles. In addition, particles at 
the maximum energy "escape" upstream from the shock system during acceleration and carry 
away energy and pressure, thus further softening t he equation of state. Th is escape is included 
implicitly in the semi-analytic DSA calculation (see iBerezhko &: Ellisonlll999l ) and is distinct from 
the diffusion of particles into and out of a particular shell. Since we use Lagrangian coordinates 
we ignore the diffusion of particles into and out of a shell after the initial acceleration. This is a 
reasonable approximation for low energy particles, but becomes less accurate as the particle energy 
and diffusion lengths increase. 

The softening of the equation of state means that compression ratios well in excess of four 
can be produ ced in non-radiative, collisionless shocks further modifying the SNR evolution (e.g., 
Eichleiill984l ). The nature of nonlinear DSA, on the other hand, uniquely predicts that a concave 
superthermal particle spectrum, hardening at higher momenta, will be produced. Here the energy 
that goes into the highest momenta particles comes from the shock-heated thermal population. 
Thus, the fundamental changes that occur with efficient DSA are larger compression ratios, lower 
temperatures of the shocked gas (compared to shocks without efficient DSA), and concave energetic 
particle spectral^ 

An important question for SNRs is whether or not particle acceleration occurs at the RS. 
While it is likely that, due to expansion, the magne tic field interior to the RS is too small to 



allow acceleration, there have been r eports of radio (jGotthelf et al 



2001 



DeLanev et al.l l2002l ) 



and even X-ray synchrotron er n ission (jRho et al.ll2002l ) at the RS in some SNRs. As discussed in 
Ellison. Decourchelle. &: Balletl ()2005l ). confirmation of particle acceleration at the RS will indicate 
strong magnetic field amplification and will have far-reaching consequences. For now, however, 
we assume there is no appreciable acceleration at the RS, but note that the large changes in the 



^The concave curvature in the superthermal portions of the spectra in the middle panel of Figure [5] is 
present but small. The curvatu re becomes more pronounced with greater acceleration effici ency (see Fig . 1 



Ellison. Berezhko. fc Barinj 200C ) and m ay actually have been detected in SNR spectra (e.g., Reynolds fc Ellison 



19921 : 



Jones et al. 



20031 : lAllen et allboosl) 
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Fig. 2. — The top and middle panels show radial distributions of the plasma temperature and 
density for SNRs with different shock acceleration efficiencies at the FS. The bottom panel shows 
p^V where each horizontal line indicates the value summed over the shocked ejecta or shocked 
CSM regions. In all panels, the heavy-weight solid curves (model Dl in table) are for efficient DSA 
(crei — 63%), the dashed curves (model B2) are for moderate DSA (erei — 36%), and the dotted 
curves (model Al) are for test-particle (erei 0). In the top panel, the light-weight solid curve 
(model D2) is an efficient acceleration case (crei — 63%) where £^SN has been chosen to produce the 
same FS radius as the TP example with £^SN = 10^^ erg In all cases, no DSA occurs at the RS. 
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density and temperature that occur with efficient DSA may make it possible to use observations of 
thermal X-rays to set limits on the amount of acceleration. 

We use as inputs to our thermal models a set of hydrodynamic models in which we vary 
the ambient magnetic field Bq, the ambient density uh, and the acceleration efficiency erei- The 
acceleration efficiency, erei, is defined as the percentage of energy flux crossing the shock (in the 
shock rest frame) that ends up in relativistic particles, and it is uniquely determined by Xmy The 
acceleration efficiency includes particles in the particle distributions downstream from the shock as 
well as the high-energy particles that escape upstream from the shock. All models are i nitialized in 



a sirn ilar way typical of Type la supernovae: we use an exponential density profile (e.g.. iDwarkadas 



2000) for the SNR ejecta, an ejecta mass Mej = 1.4 Mq, and an ejecta kinetic energy E'sn = 10^^ 
erg. We assume that the ambient magnetic field and density are uniform and evolve the models for 
a time isNRi typically taken to be 500 

In Figure [2] we show the plasma temperature, density, and volume emission measure (i.e., 
p^V) with varying injection efficiencies for DSA. In all cases, except for the light-weight curve 
in the top panel, we take nn = 0.1 cm~^, Bq = 15 fiG, tsNR = 500 yr and E'sN = 10^^ erg. For the 
light-weight curve, we use Esn = 1-4 x 10^^ erg with the same values for nu, Bq, and tsNR- Other 
parameters are listed in Table 1 as models Al, B2, Dl, and D2. 

The three sets of heavy-weight solid curves have strong DSA with ei-oi ~ 63% at isNR = 500 yr. 
The dashed curves are for a moderate acceleration efficiency with e^-d ~ 36% and the dotted curves 
are the test-particle limit with negligible energy flux in relativistic particles (erei ~ 0). These 
results clearly demonstrate the two most important effects efficient CR production has on the 
remnant dynamics: the plasma temperature behind the FS (top panel) drops by nearly a factor of 
10 between the efficient acceleration and test-particle examples, and the shock compression (middle 
panel) immediately behind the FS increases by a factor 7.7/4 ~ 20 

The quantity ^ p^V in the bottom panel (where p is the shocked plasma density and V is the 
volume of the shocked plasma) , is similar to the emission measure (EM) ; the thermal X-ray emission 
is proportional to ^ p'^V. The light-weight bar is ^ p'^V for the region between the reverse shock 
and the contact discontinuity, and the heavy-weight bar is for the region between the forward shock 
and the contact discontinuity. 

An important consequence of efficient DSA is that it results in a smaller forward shock radius at 



''A remnant age of tsNR = 500 yr is chosen mainly for computational convenience although it allows comparison 
with some young historical SNRs. 

^The compression ratios indicated in Table 1 are low compared to what is possible in nonlinear DSA (see, for 



example, lEllison et al 



20051 . 



for extreme cases with rto 



100). 



The modest values for rtot obtained here result 
from a comb i nation of restricting the acceleration efficiency and assuming laxge upstream fields. As described in 



Ellison et al. I (|2005f ). in the approximations used in our CR-hydro model, large upstream fields result in heating of 



the shock precursor and reduce the overall acceleration efficiency and rtot- 
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10^° r 11^=0.1 cm"^ 6(5=15 AiG-i 




10^ E \ , \ , \ , ; 

2 4 6 8 



radius [pc] 

Fig. 3. — Temperature, density, and ^ p'^V as in Figure [2l The various sets of curves were 
calculated at tsNR = 200, 500, and 1000 yr as indicated. In all examples, erei ~ 50% at isNR- The 
dashed curves are model A3 in Table 1. 
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a given age as compared to the test-particle limit. A comparison of the FS positions of the heavy- 
weight solid and dotted curves in Figure [2] shows that the FS radius drops by a factor of 4.6/5 ~ 0.9 
between the test-particle and highly efficient acceleration c ases with a corresp onding drop in the 
FS/CD ratio. Such an effect is observed in Tycho's SNR (jWarren et al.l l2005l ) . providing strong 
evidence for efficient particle acceleration in this young remnant. The light-weight solid curve in the 
top panel has strong DSA (erei ~ 63%), but the supernova explosion energy has been increased from 
]^g5i gj,g (heavy-weight curves) to 1.4 x 10^^ erg to obtain the same FS radius as the test-particle 
case. If efficient DSA occurs and is ignored, systematic errors in estimates for E'sN will occur. 

In Figure [3] we show the shocked plasma temperature, density, and ^snr = 200, 500, 

and 1000 yr. These are efficient forward shock acceleration cases with ej-d ~ 50% at isNR- At 
^SNR = 200 yr, is slightly higher at the RS, but by 1000 yr, X^p^F is almost an order of 

magnitude higher in the FS. This occurs, of course, because the FS continues to sweep up CSM 
material as the remnant expands. 



Figure m shows a similar set of curves, also for erei ~ 50%, only now the age is fixed at 500 yr 
and nn is varied as shown. As expected, ^ p^V at the forward shock scales roughly as n^. At the 
RS, ^is° shows a correlation with nn but it is considerably weaker than n^. 

In addition to modifying the evolution and the temperature of the shocked gas, changes in the 
compression of the fluid should result in changes in the compression of the magnetic field, implying 
that the morphology and intensity of synchrotron emission from relativistic electrons will vary 
strongly with the efficiency of DSA and with the orientation and st rength of the magnetic field | § The 
magnetic field in the CR-hydro model is exactly as described in lEllison &: Cassam-Chenaj (|2005l ) 
and we emphasize that the DSA model used here does not explicitly include the magnetic field. 
What this means is that instead of including the magnetic field self-consistently in determining the 
shock structure, the far upstream fi eld strength, Bn, is used to calculate a rate of Alfven heating in 
the shock precursor (as described in lBerezhko fc Ellisonlll999l ) which in turn modifies the subshock 
compression and overall acceleration efficiency. An ad hoc field comp ression model is assumed to 
determine the downstream field (e.g.. lEllison &: Cassam-Chenaill2005l ) and this compressed field is 
used to calculate synchrotron emission. 

Significantly, our model does not include magnetic field amplification even though this is now 



believed to be an important effect in SNR shocks (e.g., IVink &: Laming 



2003 



V51k et al.ll2005al ). 



Field amplification, as opposed to simple compression, comes about when backstreaming energetic 
particles transfer en ergy to magnetic turbulence to create a nonlinear situation where IS.B/B S> 1 
(|Ben fc Lucelj 120011). Efforts are n o w underway to incorporate this pro c ess into nonlinear models 



of D SA (i.e., lAmato &: Blasi 



2006 



Vladimir ov. Ellison. &: Bykov 



2006 



Blasi. Amato. &: Caprioli 



20061 ). We do not attempt to model such amplification here, but (in a fashion similar to, for 



^We note that the ori e ntatio n of the magnetic field is not included explicitly in our model. As discussed in 



Ellison fc Cassam-Chenail (|2005l 'l. for the purposes of calculating the evolution of the field, we assume the shocked 



field is fully turbulent. 
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radius [pc] 



Fig. 4. — Temperature, density, and p'^V as in Figure [2l The various sets of curves were 
calculated with different nn as indicated. In all ~ 50%. The dashed curves are model 

A3 in Table 1. 
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example. I Volk et al.ll2005al ). we investigate the effects of the associated strong postshock magnetic 
fields by considering the case with a strong (60 /iG) upstream field which, upon strong compression 
under conditions of efficient particle acceleration, in some ways mimics the amplification process. 

In Figure [5] we show electron and proton spectra for two cases, one with test-particle 
acceleration, (erei ~ 0; dashed curves, top panel) and one with efficient DSA, (erei ~ 63%; solid 
curves, middle panel). For this example, the specific parameters are not important other than to 
emphasize that the only difference in input between the two models is the injection efficiency at 
the forward shock Xin]- The spectra shown are integrated over the region between the CD and 
FS. Both the spectral curvature and the shift of the thermal peak to lower momenta with efficient 
DSA are evident. In the bottom panel of Figure [5] are the synchrotron and bremsstrahlung spectra 
generated by the electron distributions. The two vertical lines show the 1-10 keV range relevant 
for X-ray observations. The important points to be made here are: 

(i) The spectral curvature in the synchrotron emission is far too small to show up in the limited 
X-ray range, although it is important for i natching radio to X-ray fiuxes since nonlinear effects 



incre ase the X-ray/radio ratio significantly (jBaring et al.lll999l : Ellison et al.ll200ll : iBerezhko et al 
2OO2I ): 



(ii) In the models shown, the synchrotron spectral shape in the critical 1-10 keV range is 
determined by the shape of the turnover in the electron distribution near the maximum momentum. 
The momentum and shape of this turnover is determined by Ocut radiation losses (synchrotron and 
inverse-Compton losses off the primordial background radiation) and adiabatic expansion as well as 
the maximum momentum, Pmaxj the shock can produce. Both the rate of radiation losses and Pmax 
depend on the momentum dependence of the diffusion coefficient, k, near Pmax, which is unknown. 
To our knowledge, all serious models of S NRs currently in use, including ours, assume that k ex. p at 



200e 



Pmax, but this need not be the case (see lAmato Sz Blasil l2006l : I Vladimir ov. Ellison, fc Bykov 
Blasi. Amato. Sz Capriolil 120061 ). If k near Pmax has a momentum dependence other than oc p, the 



shape of the synchrotron emission near 1 keV can be substantially modified. 

(iii) The bremsstrahlung continua without emission lines in the bottom panel (light-weight 
curves) illustrate the effects on the thermal emission. First, note that the particle spectra in 
Figure [5] are "complete" in that they contain all of the mass in the shocked gas. For this example, 
we assume the ambient plasma is fully ionized hydrogen and helium with 10% helium by number and 
that = = T2a, that is, the shocked plasma of electrons, protons, and a-particles comes into 
equilibrium immediately behind the shock. In this example, the thermal continuum is dominated 
by the synchrotron [(heavy curve)/(light curve)] for both the TP and efficient cases. However, the 
synch/brems ratio at 1 keV drops from 10'^ in the efficient case [(heavy solid curve) /(light solid 



^The parameter Ocut is described in detail in lEllison et alj (2004) (see eq. 2 in that paper). Simply stated, when 
acut = 1, an exponential turnover is imposed on the particle spectra near pmax. This turnover can be important for 
modeling synchrotron X-rays but the actual shape of the turnover can't be determined with the semi-analytic code 
used here. 
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Log^Q Energy [keV] 

Fig. 5. — The top two panels show electron and proton spectra, p'^f{p), for TP DSA (erci ~ 1%, 
dashed curves, ~ model Al) and efficient DSA (erci ~ 63%, solid curves, model Dl). The bottom 
panel shows synchrotron and bremsstrahlung spectra for these two cases: solid curves - efficient 
DSA, dashed curves - TP. In all cases, Bq = 15 /xG and nu = 0.1 cm~^. There is a large variation 
in the synch/brem ratio in the X-ray range between these two cases. The up-down arrow in the 
top panel indicates that the normalization of the power law, relative to the thermal distribution, 
is arbitrary in the test-particle case other than that the power law is low enough to contain an 
insignificant fraction of the total energy. The up-down arrow in the middle panel shows the arbitrary 
parameter (e/p)^! = 0.01, typical of the ratio observed in galactic cosmic rays. The thermal electron 
distributions shown assume instant temperature equilibration with nrotons behind the shock. 
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curve)] to ^ 50 in the TP case [(heavy dashed curve)/(hght dashed curve)]. The shift to a lower 
electron temperature with efficient DSA is also evident in the bremsstrahlung curves. 

Actual SNRs often show a mixture of thermal and nonthermal X-ray emission (Kepler's and 
Tyc ho's SNRs show a mixture, whereas in G347.3— 0.5, the X-ray emission is wholly nonthermal. 



Slane et al.lll999l ). We now describe how we use the above hydrodynamical models to calculate 



e.g. 

the thermal X-ray emission which is expected from cosmic-ray modified shocks. 



2.2. NEI Model of Thermal Emission 



To properly predict the thermal X-ray emission produced by a shock-heated gas, it is necessary 
to determine the ionization state for each element in the gas. Particularly for the low densities 
found in SNRs, this ionization state can differ significantly from that expected for a gas in 
collisional ionization equilibrium at the associate d electron temperature of the postshock gas, Tg 



(Hamilton et al 



1983 



Hamilton &: Sarazinlll984l ). The approach to an equilibrium state for the 



ionization depends upon both the electron density ng and the time over which the gas has been 
ionizing (i.e. the time since it was shocked). The non-equilibrium i onization (NEI) state for such 
an evolving plasma has been calculated by a number of authors (e.g.. iHughes &: Helfandlll985l : lltoh 



1979 



Borkowski et al. 



2001 



the latter of whom have implemented the calculation in a routine in 
the XSPEC software package. Here the ionization state is calculated for given values of the plasma 
temperature and so-called ionization timescale r = rigt, and the ionization values are then coupled 
to a plasma emissivity code to produce the expected spectrum. 

To investigate the X-ray spectra from the shock-heated gas in an SNR under conditions of 
efficient particle acceleration, we have coupled the values for the evolution of the temperature and 
density in the postshock gas, determined from hydrodynamical simulations, with the NEI model in 
XSPEC. Specifically, for a remnant of a given age, we use the values of r and Tg, along with the 
volume emission measure nenfjdV/ {An D'^) for each radial shell from the simulation as input to the 
NEI model {D is the distance to the remnant). We consider two cases for the heating of the electrons 
in the postshock re gion: instant e quilibration with the ion temperature, and equilibration through 
Coulomb collisions (ISpitzerlll968l ). For both cases we include the effects of adiabatic expansion on 
the temperature components. We note that our line emissivity calculations assume a Maxwellian 
distribution for the exciting electrons and ignore the possible excitation by superthermal electrons 
or protons. When particle acceleration is efficient, the low energy particles may, in fact, have a 
distribution that differs somewhat from a Maxwellian although there is no reason this difference 
should be significantly more than in the case with inefficient acceleration. More importantly, 
superthermal particles may contribute to th e ionization, part i cularly for the iron K-g line. These 
effects are just starting to be considered (e.g.. lVink et al.lll997l : lPorquet et al.ll200ll ) and are beyond 
the scope of this paper. 



The ionization rate for a given species is highly dependent upon the electron temperature. As 
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Fig. 6. — The top panel shows the electron and proton temperature evolution for various shells in 
the region between the CD and FS assuming heating by Coulomb collisions. The left end-point of 
each line indicates the SNR age at which the shell was formed. The lower two panels show the total 
K- and L-shell ionization rates for several X-ray emitting ions. Ions such as Mg-K and Ar-K are 
not shown but have similar profiles. In panels a and c, nn = 1 cm~^, Bq = 15 /uG, and erci — 36% 
(model B3, Table 1). In panel b, nn = 0.1 cm~'^, Bq = 15 fiG, and ei-ci — 36% (model B2). 
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shown in the top panel of Figure [U both the electron and ion temperatures evolve as the SNR 
ages, due to Coulomb heating of the electrons and adiabatic expansion (cooling) of the mass shell. 
Therefore, the ionization state at a particular time for any given shell is dependent upon the 
ionization history in the shell, and is not a simple function of the current temperature and summed 
r. While we defer a complete treatment of the ionization directly in the hydrodynamical simulations 
to a future paper, we find that there is a suitable parameter space over which the ionization rates 
and temperatures are roughly constant in time (except for when the shell is initially shock heated) 
and, thus, for which the ionization rates for X-ray emitting ions are roughly constant (panel b in 
Figure E]). In such a case, the true ionization age of a mass shell is very close to ngtage) since the 
ionization versus time is roughly constant. In the cases where the ionization rates are not nearly 
constant with time (c.f., panel c in Figure[6]), the ionization rates for these states (i.e., Ni-K, Fe-K) 
are well below the ionization rates for the dominant line emitters (e.g., Fe-L). We note that in the 
case of shocked ejecta, this approximation will no longer be valid, since there will be more free 
electrons per ion due to the metal rich nature of the plasma, and this number will evolve as the 
ejecta ionizes. We defer the treatment of the X-ray emission from the reverse shock to a future 
publication. 

XSPEC requires the following inputs: the electron temperature Te, the ionization age net, the 
ion density rij, and the shell volume. We assume an arbitrary distance D oi 1 kiloparsec. The ion 
density is related to the hydrodynamical gas density p by = p/{fima_^^), where fi is the mean 
molecular weight for a cosmic abundance plasma, taken here to be 0.6. We then generate a NEI 
model for each mass shell. The models are summed together to produce an integrated thermal 
spectrum, such as that shown in Figure [71 These models can then either be folded through an 
instrument response to produce a simulated spectrum, or added to the nonthermal X-ray spectrum 
to produce a complete (thermal and nonthermal) X-ray spectrum|f| 



3. RESULTS 

Calculations of SNRs including particle acceleration involve a number of parameters. For all of 
our results here, we assume that the unshocked CSM is uniform (e.g., no pre-SN wind is present), 
and define the following environmental parameters for our spherically symmetric calculation: (i) 
initial kinetic energy of the ejecta material, i?sN) (h) mass of the ejecta, Mgj, (iii) age of the SNR, 
tsNR, (iv) proton number density of the unshocked CSM, nu, (v) magnetic field strength of the 
unshocked CSM, Bq, and (vi) composition of the CSM. 

Parameters and assumptions for the CR-hydro calculation include: (i) the spatial profile for 
the ejecta material, (ii) the injection efficiency for DSA, Xinj, (hi) the assumption for the scattering 



In the future, we intend to produce XSPEC table models from these combined spectra which can then be used 
to fit X-ray observations of SNR. 



- 18 - 



mean free path, i.e., A = T]mip'''g, (iv) the cutoff parameter for particle spectra near Pmax, ctcut, 
and (v) the assumption for compression of magnetic field behind the shock. The mean free path 
and ?7mfp only enter the calculation in the determination of Pmax- The maximum momentum is 
determined when the upstream diffusion length is equal to some fraction of the shock radius or 
when the acceleration time equals the shock age, whichever gives the lowest Pmax- For all of the 
models here, the fraction of the shock radius that truncates the acceleration is 0.05. 

Parameters and assumptions for the NEI calculation include: (i) the assumption for electron 
heating behind the shock, (ii) the evolution of the plasma ionization, and (iii) the chemical 
composition of the CSM. 

In all of the examples that follow, we keep the following parameters and assumptions fixed: 
SsN = 10^^ ergjfl Mej = 1.4 M0, an exponential density profile for the ejecta mass, ?7infp = 1, 
c^cut = Ij and for the magnetic fiel d configuration ac ross the shock, we assume that the field is fully 



turbulent upstream and, following I Volk et al.l (|2002l ). set the immediate downstream magnetic field 



^2 = VV3 + 2r4j3i?o , (1) 
where rtot is the shock compression ratio0 

An example product is given in Figure \7\ where we show emission in the X-ray energy range 
for a particular set of parameters. For this illustration we have nn = 0.1 cm~'^, Bq = 15 /.tG, and 
the results are shown at tsNR = 500 yr. Here and elsewhere, we only calculate line emission from 
the forward shock, summed over the region between the FS and contact discontinuity, and do not 
adjust the spectra for any type of instrument response. We assume a solar composition for the 
CSM and, most importantly, we assume that the ionization age r of the shocked CSM is given 
solely by the product of the final electron density and the age of the mass shell. 

In the top panel of Figure [7] we show a test-particle example where no significant population of 
relativistic particles is produced. The middle and bottom panels show the result for efficient DSA, 
where Xinj has been chosen so approximately 50% of the shock energy flux goes into relativistic 
particles, mainly protons0 The dashed curve in the middle panel is the synchrotron emission 
produced by TeV electrons. Except for a few strong lines, this dominates over the thermal emission. 
The bottom panel shows the summed thermal and nonthermal spectra from the middle panel. 
A direct comparison with the top panel shows the difference that results by changing a single 
parameter, the assumed efficiency for DSA. For this example, the shocked electrons were heated 
via Coulomb collisions with the hot, shocked protons. 



®We use EsN = 1.4 x 10*^ in a single example shown in Figure (2] 

^°T he details of this assumption and for the evolution of B behind the shock are given in Ellison fc Cassam-Chenai 



( 20051 ). 

^^In all of our examples, the efficiency is calculated at the end of the simulation. During earlier times, the 
instantaneous acceleration efficiency can be greater or less than the stated value. 
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Fig. 7. — The top panel shows thermal X-ray emission with no DSA (crei = 0, model Al). The 
middle and bottom panels show spectra when the percentage of energy placed into relativistic 
particles by the FS is ^ 50% (model A3 in Table 1). In all cases, only emission from the FS is 
shown and it is summed over the region between the FS and CD at tsNR = 500 yr. The dashed curve 
in the middle panel is the synchrotron flux from electrons accelerated by the FS and the bottom 
panel shows the total emission. For this example, Coulomb heating is assumed for electrons and 
the mean downstream temperatures are indicated. 
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Bq = 15 /^G, = 0.1 cm 
Instant Equil. Coulonnb heating 




Energy [keV] Energy [keV] 

Fig. 8. — Thermal emission for examples where the efficiency of DSA at the FS is varied as indicated. 
The top panels are model Al, the middle panels are model A2, and the bottom panels are model 
A3 in Table 1. The emission is summed for the region between the CD and FS and emission from 
the RS is not included. The dashed curve in each panel is the synchrotron emission from electrons 
accelerated at the FS. Panels on the left assume instant equilibration between electrons and protons 
and panels on the right assume Coulomb heating of electrons. All parameters for the various models 
are identical other than erei and the electron heating. 



- 21 - 



Bq = 15 fiG, Eff = 36% 
Instant Equll. 



Coulonnb heating 
' ' ' ' I 




0.1 r Hh = 1 cnn 



0.01 



b I I I I 



<Tg> ~ 8.0 keV 



_l I I I 1_± I I I I l_L 



<T.> ~ 1.1 keV I 



_i I I I i_ 



0.5 1 2 5 0.5 1 2 5 

Energy [keV] Energy [keV] 



Fig. 9. — Examples similar to Figure [8] where erci is set at 36% and nn is varied as indicated. The 
top panels are model Bl, the middle panels are model B2, and the bottom panels are model B3 in 
Table 1. 
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Fig. 10. — Examples similar to Figure [8] where Crei is set at 36%, nn = 0.1 cm"'^ and the unshocked 
magnetic field is varied as indicated. The top panels are model CI, the middle panels are model 
C2, and the bottom panels are model C3 in Table 1. 
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In Figures [SHlOl we show results for models A, B, and C in Table 1. We present the results 
for the two extremes for electron heating (instant equilibration and Coulomb heating) to stress the 
obvious differences in the thermal continuum which arise from assuming high electron temperatures. 
We also vary the acceleration efficiency, erei, the CSM density, nn, and the magnetic field, Bq. In 
Figure [HI we have varied erei between and 50% while keeping Bq = 15 //G and nn = 0.1 cm~^. 
Figure [9] shows the effects of varying the density ?ih while holding erei and Bq constant. Finally, in 
Figure fTOl we vary the magnetic field while holding erci and nn constant. In these figures we plot 
the nonthermal synchrotron emission as a dashed line and the thermal spectrum as a solid line but 
do not plot the summed spectrum. 

We chose to focus mainly on varying the ambient density, acceleration efficiency, and magnetic 
field because these parameters have a large influence on the SNR evolution. Increasing the upstream 
density results in higher volume emission measures and thus higher thermal continua. Increasing 
the acceleration efficiency changes the equation of state of the shocked CSM resulting in lower 
post shock plasma temperatures and higher densities. Finally, varying the upstream magnetic field 
strength influences the acceleration process because it is assumed in our model that energy in 
shock-accelerated particles in the precursor can be transferred into magnetic turbulence and then 
into heat. This heating of the precursor w eakens the subshock and results in less efficient particle 



acceleration (see iBerezhko fc Ellisonlll999l . for a full discussion). For the cases shown in Figure [TO 



the injection efficiency, Xinj has been adjusted to produce erei — 36% as Bq is changed. 
General results from Figures El [9l and [lO] include: 

(i) synchrotron dominates thermal emission except for the TP examples (top panels in Figure [8]) 
and the nn = 1 cm~^ examples in the lower panels in Figure [H 

(ii) In the TP case (top panels in Figure Ej), the differences between instant equilibration and 
Coulomb heating are large, with the average electron temperature for instant equilibration ~ 35 keV 
and only ~ 0.73 keV for Coulomb heating. However, once synchrotron becomes important with 
Crci ^ 10%, the signature of electron heating is difficult to discern against the synchrotron continuum 
except for high nn (i.e., the nn = lcm~^ cases in Figure [9]). 

(iii) The cases involving instantaneous equilibration produce models with anomalously high electron 
temperatures, even when DSA is efficient. While these simulations cannot rule out instantaneous 
electron heating as we are not comparing them to observations of SNR forward shocks, fits to X-ray 
spectra do not generally predict such high electron temperatures in these regions. 

In the top panel of Figure [11] we show the electron and proton spectra for the three instant 
equilibration examples of Figure [TOl The bottom panel shows the synchrotron spectra produced by 
these electrons. The three cases have erei — 36% even though the proton spectra look quite different. 
The major reason for this is that e^ei includes the energy in particles that escape upstream from 
the shock and which is not included in the plotted spectra. The spectra shown are summed over 
the interaction region and the acceleration efficiency varies during the evolution of the remnant. 
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= 0.1 cm ^, DSA Efficiency = 36% 




Log^Q keV 

Fig. 11. — The top panel shows particle spectra summed over the SNR at the end of the simulation, 
i.e., isNR = 500 yr. The bottom panel shows the synchrotron emission from these electrons. All 
models have the input same parameters except for Bq which is 3 nG for the dotted curves (model 
CI), 15 /xG for the dashed curves (model C2), and 60 /xG for the solid curves (model C3). 
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Fig. 12. — Pion-decay and inverse-Compton emission for a range of nn and Bq. In the top panel, the 
heavy-weight curves are pion-decay, the hght-'weight curves are inverse-Compton, and e^ei = 36% 
and Bq = 15 fiG in all cases. The strong dependence of pion-decay on ambient density nn is 
evident. The middle panel sho'ws IC and the bottom panel sho'ws pion-decay for nn = 0.1 cm~'^ 
■with Bq varying from 3//G (solid curves) to 15 /xG (dashed curves) to 60 //G (dotted curves). For 
comparison to the ir^, we show in the bottom panel the IC emission for Bq = 60 (light-weight 
dotted curve). The particle distributions producing the emission in the bottom two panels are 
those shown in the top panel of Figure [TTl 
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As Bq increases, the maximum proton momentum increases in approximate proportion to i?nl I 
The maximum electron energy decreases with Bq, however, because synchrotron losses become more 
severe. The break in the electron spectrum for the Bq = 60 fiG case near W^mpC (solid curves) 
results because strong radiation losses lower Pmax for shells shocked early-on. When all shells are 
added together at fsNR, a steep region spanning several decades in momentum results and this 
produces a relatively flat region in the X-ray synchrotron spectrum shown in the bottom panel. 
In the 1-10 keV X-ray energy range, the synchrotron flux varies by nearly an order of magnitude 
between Bq = 3 and 15 /uG, but the shape is almost indistinguishable. 

The GeV-TeV electrons produce inverse-Compton emission as well as synchrotron radiation, 
and the protons accelerated by the FS interact and produce GeV-TeV photons via proton-proton 
collisions and pion-decay. In the top panel of Figure [12] we compare the inverse-Compton and pion- 
decay emission for the examples shown in Figure [9] where tt-h varies and B = 15 //G. As expected, 
the overall pion-decay flux increases more rapidly with nn than the inverse-Compton flux. At 1 
TeV, the ratio (pion/IC)iTeV is 0.01 for nn = 0.01 cm~^, 0.08 for tt-h = 0.1 cm~^, and 0.65 for 
nn = lcm~^. It is noteworthy that in order to obtain (pion/IC)iTeV > 1 densities greater than 
1 cm~'^ are required. At energies below ~ 100 GeV, the pion-decay emission for nn = 1 cm~^ 
becomes stronger than the inverse-Compton as it does at energies above a few TeV. The energy 
range that will be explored by the Large Area Telescope (LAT) on GLAST (20 MeV-300 GeV) 
should allow a clear differentiation between emission processes in this case. 

In the bottom two panels of Figure [l2] we show the IC and pion-decay emission for the particle 
spectra shown in the top panel of Figure [Til The radiation losses caused by the large magnetic field 
(dotted curves: Bq = 60 /iG and B2 = 310 /zG), combined with the remnant evolution, produce 
a break in the IC emission in the GeV-TeV range. A comparison of the light- and heavy-weight 
dotted curves in the bottom panel shows that the shapes of the IC and vr'^ components are less 
distinct than with lower fields and this may make it more difficult to distinguish these components. 
When Bq is large, however, the pion-decay emission extends to much higher energy and combining 
the observations of GLAST and ground-based air-Cherenkov telescopes will help in distinguishing 
between the hadronic and leptonic scenarios. Currently, HESS spectra of SNRs roughly span the 
energy range 0.1 — 10 TeV and above 1 TeV the difference between the two dotted curves in the 
bottom panel becomes significant. 



4. DISCUSSION AND CONCLUSIONS 

The X-ray emission from SNR shocks is generally a mixture of thermal continuum and line 
emission, and nonthermal continuum. These components are connected through the particle 



^^We note that this simple scaUng is almost certain not to hol d if magnetic f ield am plification is considered. A full 
discussion is beyond the scope of this paper but preliminary work (jAmato fc Bla si 2006: iVladimirov. Ellison, fc Bykov 



20061 : Blasi. Amato. fc Capriolil 2006l l suggests that pmax will not increase in proportion to Bq. 
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acceleration process, but are generally modeled independently. Here we present the first results 
from an ongoing project to model the X-ray emission in a fully self-consistent manner with the SNR 
hydrodynamics and broad-band continuum radiation. The problem is a complicated one with a 
number of uncertain parameters both for the model and the environment. In order to simplify our 
presentation and to emphasize the important aspects of our model, we concentrate on emission from 
the outer blast wave (forward shock) and use approximate techniques to model the non-equilibrium 
ionization in the interaction region between the contact discontinuity and FS. 

The results we present are intended to show broadly how the X-ray emission depends on 
important parameters including the shock acceleration efficiency, erci, the unshocked CSM proton 
number density, njj, the magnetic field strength, Bq, and two extremes of electron heating: instant 
equilibration and Coulomb heating. In summary, our results for the forward shock show: 

(i) Synchrotron emission from relativistic electrons becomes important, and even dominant, in 
the X-ray energy range for modest acceleration efficiencies over a wide range of CSM densities 
and magnetic fields. As shown in Figures [HI [H and [TOl only TP cases or those with large tt-h ^ 
1 cm~^ show thermal line emission clearly dominant over synchrotron emission. Once e^cl ^ 10% 
synchrotron emission is likely to be important for typical SNR parameters. 

(ii) The two extremes for electron heating, instant equilibration and Coulomb heating, produce 
very different thermal spectra as expected. However, it may be difficult to distinguish these in 
cases where synchrotron radiation dominates. 

(iii) Changes in the ambient magnetic field, Bq, have a large effect on the synchrotron emission 
in the X-ray range but these effects tend to saturate in terms of overall flux for large Bq (bottom 
four panels in Figure [TO]) . The saturation is a result of the fact that large magnetic fields tend to 
dampen DSA, at least in the approximate DSA model used here. 

(iv) As has been recognized for some time, broad-band models and fits are critical for constraining 
important SNR parameters; we illustrate this in Figures [11] and [T2l In homogeneous nonlinear 
models, such as ours, changes in parameters that modify emission in one energy band will modify 
emission from radio to TeV 7-rays. Thus, if large magnetic fields produce hard synchrotron spectra 
in the X-ray band, these same electrons will produce strong radio emission (bottom panel of 
Figure [TT]) . If particular ambient densities are required to match thermal X-ray observations, 
these densities strongly influence the (pion/IC)iTeV ratio (top panel in Figure [T2]) . 

(v) As has been emphasized by many authors, direct evidence for CR ion production in SNRs can 
come from identifying pion-decay emission at GeV-TeV energies. The interpretation of GeV-TeV 
observations is complicated however because pion-decay photons from proton-proton interactions 
and IC emission from TeV electrons scattering off of the cosmic background radiation are expected 
to produce similar fluxes in the GeV-TeV energy range. It is generally assumed that the shape 
of pion-decay emission is sufficiently different from that of IC to allow a clear differentiation if a 
large enough energy range is sampled. In the bottom two panels of Figure [12] we show spectra 
for a range of Bq and it is clear that the shapes of the IC spectra in the GeV-TeV region depend 
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fairly strongly on the magnetic field with the Bq = 60 fj,G case (where strong radiation losses and 
evolution produce a break in the underlying electron spectrum) closely mimicking the pion-decay 
shape in the 1 GeV-1 TeV range. For large fields, however, significant differences in shape and 
intensity occur at both lower and higher energies highlighting the importance of air-Cherenkov 
observations above 1 TeV and for furure observations by GLAST below 1 GeV. 

The model we have presented is, we believe, the first to self-consistently determine the 
synchrotron continuum with thermal emission in the X-ray energy range. However, a number of 
steps remain before realistic models are produced that can be used to interpret X-ray observations. 
These steps include: (i) modeling of the reverse shock, which requires an accurate description of the 
ejecta composition and spatial distribution; (ii) calculating the ionization rates for material behind 
both the forward and reverse shocks, including time dependence and variations in temperature and 
electron density; (iii) line-of-sight projections; and (iv) incorporating the effects of the instrument 
response for general X-ray spectra without restriction to specific models such as those currently 
used in, for example, XSPEC. 
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